Coulomb corrections to the Delbriick scattering amplitude at low energies 
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In this article, we study the Coulomb corrections to the Delbriick scattering amplitude. We con- 
sider the limit when the energy of the photon is much less than the electron mass. The calculations 
are carried out in the coordinate representation using the exact relativistic Green function of an 
electron in a Coulomb field. The resulting relative corrections are of the order of a few percent for 
scattering on for a large charge of the nucleus. We compare the corrections with the corresponding 
ones calculated through the dispersion integral of the pair production cross section and also with the 
magnetic loop contribution to the g-factor of a bound electron. The last one is in a good agreement 
with our results but the corrections calculated through the dispersion relation are not. 



INTRODUCTION 

The elastic scattering of photons by an external 
Coulomb field (the so called Delbriick scattering [l|) is 
one of nontrivial predictions of quantum electrodynam- 
ics. In the perturbation theory, the Delbriick scattering 
amplitude begins from the second order in Za (Z\e\ is 
the charge of the nucleus, a = e 2 w 1/137 is the fine- 
structure constant, we put c = 1, % = 1). Therefore, 
significant efforts have been made to calculate the ampli- 
tude for the arbitrary scattering angles and energies even 
in the lowest-order Born approximation. The results of 
these calculations and the detailed bibliography can be 
found in the report Q. 

To calculate the Delbriick scattering amplitude for 
Z ^> 1 it is necessary to take into account Coulomb field 
exactly. The expression for the amplitude exact in Za 
has been derived in Ref. Q without any additional as- 
sumptions, but numerical results have not yet been ob- 
tained because it is fairly cumbersome. 

Considerable progress in the calculation of the 
Coulomb corrections to the lowest-order Born approxi- 
mation has been achieved for the case of the photon en- 
ergy ui much large than the electron mass m e and either 
at small scattering angles A/u> <C 1 (A = |ki— k 2 |, where 
\ii and k 2 are the momenta of the photon Jn_the_ initial 
and final states, correspondingly) Refs. 
at large momentum transfer A/m e 
It turns out, that the Coulomb corrections strongly de- 
crease the Delbriick amplitude in comparison with the 
lowest-order Born approximation (the theoretical results 
and the corresponding experimental data are reviewed 
in detail in the report 11 1 and the tables [13]). At the 
moment, the minimal photon energy at which Delbriick 
scatterin g is exp erimentally observed is u> = 889 KeV 
Refs. 13, 14, [lj| [IH, [13] but the corresponding energy 
for the Coulomb corrections is u — 2754 KeV Ref. [18f 
only, i. e. above the electron-positron pair production 
threshold. 

We aim here to calculate the Coulomb corrections at 
low energy to -C m e . These corrections have not yet 
been investigated neither experimentally nor theoreti- 




cally. Nevertheless, they are closely connected with the 
Coulomb corrections to the pair production cross section 
due to the dispersion relation Refs. lj|, 2(| and, as we will 
show below, with the magnetic-loop contribution to the 
g-factor of a bound electron Ref. [2l| . We think that the 
Delbriick scattering amplitude at low energy is useful to 
estimate the Coulomb corrections for both phenomena. 

The amplitude is calculated in the coordinate represen- 
tation with the help of the Green function of an electron 
moving in a Coulomb field. The structure of this paper 
is the following: in Section 2, we provide all necessary 
information about the Green function and the general 
parametrization of the Delbriick scattering amplitude. In 
Section 3, we show that the calculation in the coordinate 
representation reproduces the result of the lowest-order 
Born approximation derived in the momentum represen- 
tation. We point out the difficulties specific for a cal- 
culation in the coordinate representation also occurring 
during the calculation of the Coulomb corrections. The 
results for the Coulomb corrections are given in Section 
4. We also provide the simple parametrization of their 
dependence on Z. In Section 5 we compare our results 
with those obtained via the dispersion relation. The es- 
timated value for the magnetic-loop contribution to the 
(/-factor of a bound electron is given in Section 6. 



DELBRUCK SCATTERING AMPLITUDE 

We parameterize the Delbriick scattering amplitude as 
follows: 



A = eWe*W IP" ( w ,k 1; k 2 ,Z). 



(1) 



IF" (w,ki,k 2 ,Z) = 



a {Zo>Y 



+h (w, ki, k 2 , Z) [tu 2 g^ - to (n"A£ + n"*£) 

+ n' i n v k 1 -h 2 ]} (2) 



2 



where k% — (w,ki),fe2 = (w,k 2 ) are the 4-momenta of 
the photon in the initial and final states, correspond- 
ingly, e^ 1 * 2 ) are the polarization vectors, the 4-vector n 
is defined as ki ■ n — k 2 ■ n = uj, fi and f 2 are the form 
factors to be calculated. The main purpose of this article 
is to calculate lim| k |^ /i,2 (0, k, k, Z). 

In a point- like charge approximation (Coulomb field), 
the polarization tensor IP" has the following form: 

IF" (. . . , Z) = IP" (. . . , Z) - IP" (. . . , 0) , (3) 

n A " y (a),ki,k 2 ,2') = ict J d 3 r 1 d 3 r 2 exp(ikiri - ik 2 r 2 ) 

x f ^Sp{ 7 ^G(r 1 ,r 2 |6) 7 "G(r 2 ,r 1 | £ -c)}, (4) 
J c 



where G (ri, r 2 |e) is the Green function of an electron in 
a Coulomb field. The contour of integration over e in the 
expression (0| goes from — oo to oo so that it is below the 
real axis in the left half-plane and above the real axis in 
the right half-plane. It is convenient to turn the contour 
along the imaginary axis. In this case the Green function 
takes the following form (see Ref. [22|): 



^ 00 />oo 

G(ri,r 2 I ie) = } / ds exp(2iZa s (e/p) — p(r± + r 2 ) coths) 

Airr x r 2 p f^Jo 

X { [#+! I' 2v {y)Bi + R-U 2v (y)Ai (70 ie + m) + Zaj° [im (n x + n 2 ) + p R+ coths] I 2v (y) B, 

1 hu(y)\, 



where 



p 2 * l2 2 ( fi i + "2) Bi + pcoths(n x - n 2 ) lA t 
2 sinn s 



R± = l±mn 2 ±i£ (ni x n 2 ) , E fc = ie ijk [7*, -f j ] /4, n (1>2) = -yn (1>2) , 



dx 



Bi = — (Pi(x) - Pi-i(x)) , cc = n 1 n 2 , 
dx 



v = Ji 2 ~{z a y 



y = 2p^Jr\r 2 j sinhs, p — y m 2 + e 2 . 



(5) 



(6) 



For the sake of convenience, we calculate the time-time 
component of the polarization tensor and the trace of the 
spatial components separately: 

(k, Z) = (n„n„ - gixv ) IP" (0, k, k, Z) , (7) 



I!* 00 ) (k, Z) = n^W" (k, Z) . 



(8) 



Substituting the parametrization of IP" ^ in the right- 
hand side of the Eqs. (7J) and JSj) yields the relation be- 
tween !p J ),]P 0CI ) and the form factors /(i,2) : 



n («) (k! z) = 2^^k 2 /i(o,k,k, 

mi 



Z) (9) 



n<°°> (k,z) = -^!k 2 (/ l( o,k,k, 



z) 



+/ 2 (0,k,k, Z)). 



(10) 





(a) (b) 
FIG. 1: lowest-order Born approximation 

LOWEST-ORDER BORN APPROXIMATION 

The diagrams of the second order of the perturbation 
theory in Za are depicted in Fig.[T] Their contribution 
were calculated in Refs. [23, 24 1. We aim here to demon- 
strate that the calculation of this diagrams in the coordi- 
nate representation reproduces the result of the lowest- 
order Born approximation. 

To calculate the contributions of the diagrams, see 
Fig.[TJ we expand the exponent function in the expres- 
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sion ([4]) up to the second order in |ki| = |k 2 | = |k|: of the n-th order in Za: 

2n (a)+ n ( b )-^/ Mr^-r^J^ 
x Sp{2 7 ^G(°) (ri,r 2 |ze) 7 l 'G , ( 2) (r 2 ,n|ie) 
+7 M G (1) (r u r 2 \ie)YG^ (r a ,n|ie)} , (11) 

where G^ is the contribution to the Green function ([5|) 



G (o) = _ 



( 7o ie + m) + i 



ri - r 2 



d 



Iri - r 2 | d Iri - r 2 | 



exp(- |ri -r 2 |) 
47r|ri-r 2 | 



(12) 



where 



G« = - 



Za 1 - t 2 ( ep 
1 2i-W 



«P /- ~ % 9 ip „ . <9 
7 ze + m+ -3- (ri - r 2 ) — - + — (m - n 2 ) — 
p^ 000 2 dp 

-p 



167r rir 2 [_ p 2 

1 n / " ™> , . d \ f 1 

— (*-- 7 (ftl + ft2) ^) U e e 



G ( 2 ) _ (^«) /" rfs exp(ip(ri + r 2 ) coths) < 



47rr 

. p(n - r 2 ) 



2/9 



R+B l j^-+ R-lAi ) (70 te/p + m/p) 



, (ni + n 2 ) - i coth s (ni - n 2 ) ZA; 
2 sinh s 



41 dv 



es 

+ 5 



— T + 7 — (ni + n 2 ) - coth s 
2p 



hu{y) 



h [y cos - 



e-" CT log + eP° T (0, p (1 + a)) e^ a Y (0, p (1 - <r)) 

1 — a 



T : 



2j/ 2 (70 ie/p + m/p) - iy 2 coths (nj — n 2 ) — « 



p(n - r 2 ) (ni + n 2 ) ?yd 
sinh 2 s cos 2 § <9y 



I [y cos ■ 



(13) 
(14) 
(15) 
(16) 

(17) 
(18) 



where T (a, 2) = i a_1 exp (— i) eft is the incomplete 
gamma function. The following notations have been also 
introduced in the expressions (|12H18p : 

P = p(n+r 2 ), „ = \E±-lA, (19) 

ri + r 2 

, n-r 2 4> /l + nin 2 \ 1/2 

t= — , cos o= o ' 20 

n + r 2 2 V 2 J 

The contribution to the Green function G^ 2 ' of the second 
order is split into two parts: the contribution correspond- 
ing to the expansion of a Bessel function over an index 
(the part of the expressions (fT4|l , (fl~5|) separated by the 
square brackets) and the contribution in which we sub- 
stitute the indices of the Bessel functions for 21 and sum 
over I explicitly. The latter can be called conditionally 
quasiclassical contribution because it corresponds to the 
contribution of the large angular momenta I 3> Za. 

There are some points to be made. First of all, the con- 
tribution of each diagram Fig.[T]is infrared divergent, i. e. 
it is divergent at large distances. Since the divergence is 



canceled between the diagrams Fig.flja) and Fig.QJb), 
the contribution each of them depends on the regular- 
ization of this divergence. However, the contribution of 
separated terms (for example, the quasiclassical contri- 
bution or the contribution proportional to d v l 2v (|14[) . 
p5p) and even the presence of the divergence depends 
on the order of the iterated integration over the spatial 
variables r , and the inner variables of the Green functions 
- "proper times" Sj. An iterated improper integral, the 
value of which depends on the order of integration, we 
call a conditionally convergent iterated integral. The ex- 
ample of such integral, that arises during the calculation 
of the diagram Fig.fljb), is given in Appendix. To avoid 
the complications due to an explicit regularization and 
difficulties during the calculation of separated terms, one 
should fix the order of integration for all the diagrams 
and sum the contribution of each one before the inte- 
gration with respect to the last variable. We chose the 
variable t defined in (|2"D|) as the last integration variable. 
In this case, the contributions of the diagrams (Fig.[T|) 
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FIG. 2: The relative Coulomb corrections to the trace of the 
time-time component of the polarization tensor. The dashed 
curve corresponds to the fit a(Za) 2 , the solid line corresponds 
to a(Zaf + b(Za) 4 
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FIG. 3: The relative Coulomb corrections to the trace of the 
spatial components of the polarization tensor. The dashed 
curve corresponds to the fit a(Za) 2 , the solid line corresponds 
to a(Zaf + b(Za) 4 



have the following form: 
n («) = a{Za) 2 k2 



n 



_ a{Zaf 

(b) - _a 



k 2 



5 


5 




2304 


~ 128 


Jo 1 


( 19 




r?) 


VH52 







(21) 



(22) 



n (») _ 2n («) + n (ii) - a( - Za) k 2 — (23) 

The time-time component is derived in a similar manner: 



(oo) = a{Za) 2 2 59 
Born ml 2304' 



(24) 



Substituting the expressions ([S3]) and ((Ml) in © and 
yields the form factors ((31) in the lowest-order Born ap- 
proximation 

/lB= 16^72' /2B= "32?72- (25) 

As noticed above, they coincide with the results derived 
in (Hill. 



COULOMB CORRECTIONS 

Analytical derivation of the Coulomb corrections to the 
lowest-order Born approximation (|25[) is a rather com- 
plicated problem. We have calculated these corrections 
mostly numerically. To increase the accuracy of the nu- 
merical calculations we have subtracted the lowest-order 
Born approximation Eq. (|lip from the general expression 
Eq. ([3]) before any transformations. Let us denote the an- 
gular momenta and the proper times that appear in the 



Green functions in Eq. ^ as l±, I2 and si, S2. After sub- 
traction, it is convenient to change variables as follows: 
(?T,r2) — > (77 = y/rxT2, f = ri/r<i). Then we integrate 
analytically over e, x = (rir 2 )/|r 2 ||ri|, t' one by one. Af- 
ter that we also perform the analytical summation over 
li. Thus, the expression (g]) can be reduced to the sum 
over I2 and the iterated integral over s%, s 2 and 77. The 
explicit expression for the integrand is omitted here as 
bulky. Further analytical integration is only possible for 
separate terms. The term previously referred to as the 
quasiclas steal (which contains I21 instead of I 2i ,) can be 
represented as an one-fold integral or as an infinite series 
over Za. For example, the corresponding contribution to 
I!* 00 ) is the following : 



n 



(00) 



n 



(00) 



a(Za) 2 



quasicl. quasicl.Born 



/287«*-39 a _ ^(48^-15) 4 
V 18432 v ' 69120 ^ ; 

,^(15^-63) «_ ^(21^-10) „ 



3096576 



5160960 



^ (832907T 2 - 46431) 10 
245248819200 1 a) 



(26) 



This contribution is finite, i. e. the infrared divergency is 
absent in Eq. (|26)) . Nevertheless, the residual part must 
be integrated in the same order as that used to derive 
Eq. (|2"6")l . We check explicitly that the contribution con- 
taining the subtraction from the single Bessel function, 
i. e. which is proportional to 



+ 2v ~ Hi — 21 O v l2v 



= 0{Z 4 a 4 ), 



(27) 



is a conditionally convergent iterated integral. 

The results of our numerical calculations are presented 
in Figs.dJ[3]and Tab.|H 
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TABLE I: Relative Coulomb corrections 



z 


( n coo) _ ncoor 


\ /n (oo) 

1 / 11 Born 




/tt(w) 

/ U Born 


50 


1.82 • 10" 


2 


8.20 • 10" 


-3 


60 


2.69 ■ 10" 


2 


1.15 ■ 10" 


-2 


70 


3.78 • 10" 


2 


1.60 ■ 10" 


-2 


80 


5.15 • 10" 


2 


2.19 ■ 10" 


-2 



It should be noted that the quasiclassical contribution 
is of the order of 60% of the lowest-order Born approxi- 
mation for Z = 80. The residual part containing the sub- 
straction from Bessel functions reduces the corrections up 
to a few percents. This cancellation adversely affects the 
accuracy of the calculation. The authors carried out the 
calculations of the corrections independently. The accu- 
racy of the results (Tab.|I|) is better than one percent. 

Now we consider the dependence of the Coulomb cor- 
rections on Z. The first two terms in Eq. (f26|) domi- 
nate and give 97— 96% of the quasiclassical contribution 
in spite of the fact that we calculate complete series in 
Za. It comes as a surprise that the complete results can 
be adequately fitted by a biquadratic polynomial in Za 
without a free term: 



n( oo) _ n (oo) 2304 



n 



(00) 
Born 



59 



(3.22 ±0.01) • 10~ 3 (Za) 



+ (1.90 ± 0.02) • \Q- A {Zay , (28) 



PAIR PRODUCTION VIA DISPERSION 
RELATION 

In the papers (2(j, Gluckstern and Rohrlich have 
derived the relation between the pair production cross 
section in a Coulomb field and the Delbriick amplitude 
averaged over the polarizations: 



H = ^L ^-^-H0 rf "- 



(32) 



The amplitude (p} averaged over the polarizations of the 
photon has the form: 

A= 1 ( 5 iJ = -^/V. (33) 

One can find the relation between the Coulomb correc- 
tions to the form factor fi and the pair production cross 
section in a Coulomb field by using the expressions 
(we put m e = 1 in this section): 



1 



2n 2 a (Za) Ji 



a (J) 



duj'. 



(34) 



Let us check the formula ()34|1 in the Born approximation. 
Substituting Z = 82 (lead) yields: 



1 



2n 2 a (Za) f2B Ji 



-o B (<*/) 



dw' = 1 + 4 • 10" 



(35) 



n 



(«) 

Born 



n 



Born 



576 r 

^ [(6.69 ±0.17)- 10- 4 (Za) 
(3.18 ±0.54) • lO' 4 (Za) 4 



(29) 



The results of the fitting with a quadratic function 
a(Za) 2 and also the functions in Eqs. ([28)1 . (j29|) are 
shown in Figs.[2l[3] One further comment is in order. 
The coefficients at (Za) 2 in Eqs. (J2H1) and ^ have a 
magnitude one or two orders less than those at (Za)° 
in the lowest-order Born approximation Eqs. (|2"3"|) , 
If one assumes the same hierarchy between the coeffi- 
cients at (Za) 2 and (Za) 4 , then the coefficients of (Za) 4 
could not be distinguished from zero with our accuracy. 
In this case, the maximal difference between the dashed 
and solid curves in Figs.[2j[3] shows the actual accuracy 
of our calculations. Substituting Eqs. (j2"5)) , (|2T)|) in the 
relations Eqs. ([9]), (p~0|) . we obtain the form factors /(i 2) : 



fx = 



7 



1152 



3.35 • lQ- 4 (Za) 2 + 1.6 • 10~ 4 (Za) 4 , (30) 



73 
2304 



3.36- 10~ 3 (Zaf 



2.1 • 10- 3 (Za) 4 . (31) 



where ctb is replaced by the asymptotical formulae de- 
rived by Maximon in Ref. [25] for uj < 2.1: 



o B (u) = a(Za) — 

3 \ LU 



29 
960 



/ e 23 9 11 

x l H 1 e -i 

V 2 40 60 

2cj-4 
6 ~ 2 + w + 2(2w) 1 /2' 

for uj > 2.1: 
a B M = a (Za) jylog2w- — 



e 4 + 0(e 5 ) 



(36) 
(37) 



6 log 2uj - - + - log 6 2uj - log 2 2uj 



-^lo g 2w + ^ + 2C (3) 
3 6 



29 log 2uj 



9 • 256 




— log 2uj + - 
16 B 8 



27-512 



(38) 



Now we discuss the Coulomb corrections to the form fac- 
tors. Using the Eq. (|3"Tj) we have obtained the relative 
correction to the form factor in the Born approximation 



6 



o.oi 



a 0.005 
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FIG. 4: Coulomb corrections to the pair production cross 
section (Z = 82) 



(here and below all the calculations are carried out for 
Z = 82): 



h - hi 



f: 



4.9 • 10" 



2 LI 



(39) 



However, if we use the Coulomb corrections to the pair 
production cross section ac (uj) = a (uj) — as (uj) derived 
in Ref. 26] for the photon energy ui < 10 and the inter- 
polation equation derived in Ref. [27} for uj > 10, then 
the relative correction to the form factor in the Born ap- 
proximation is 



27r 2 a {ZaY f 2 B h 



a c (uj) 



duj = 2.7 ■ 10" 



(40) 



01 



This result is 20 times less than that in Eq. (f3"9"|) . The 
integrand (|4"0|) as a function of uj is shown in Fig.UJ It 
varies mainly in the region 2 < u> < 30 but there is a long 
negative "tail" for uj — > oo. The total integral is a result 
of the almost complete cancelation between the positive 
contribution for u) < 10 and the negative one for uj > 10. 
The following ratio shows it clearly: 



J 2 °° <jt (w) /uj 2 dui 
/ 2 °°|cr T (w) /uj 2 \ dio 



3.9 • 10" 



(41) 



For the integral (|4"0|) to be calculated with the sufficient 
accuracy it is necessary to derive the Coulomb correc- 
tions to the cross section with an accuracy better than 
a few percents. It is quite possible that this cancelation 
causes the discrepancy due to the lack of precision in 
the calculations of the positive part of the integrand in 
Ref. [3. 

The cause of the discrepancy could also be the inter- 
polation equation derived in Ref. [27| (the region 10 < 
uj < 30). Another interpolation formula for the Coulomb 
corrections to the pair production process is derived in 
Refs. [2^, up to terms which are of the order of l/u> 



TABLE II: Integration of the l/a;-corrections over the nega- 
tive "tail" [see Fig.Hand Eq. 



Contribution in ac when uj — » 00 — 


r, °c(uj) 


0(1) 


-0.184 


0(l)+0(l/u) 


0.068 


0(1) + 0(l/w) + 0(1/uj 2 ) 


-0.062 



and {1 / uj 2 ) log uj/ 2: 



9 uj 



a c {uj > 2) = a (aZ) 
-MZafMZa)] + Al 0g ||, 



Im g(Za) 



(42) 



where the functions /, g and fx are derived analyti- 
cally but the coefficient b is obtained by an interpola- 
tion procedure from the experimental data for ui > 30 
Refs. (3^, 31 1 . The absolute value of the approxima- 



tion formula (|4"2")) is always less than the correspond- 
ing corrections of Ref. (27| for uj > 25. The expression 
flUJ is zero when u) = 8.95 (see Fig.fJ] the correspond- 
ing value for the approximation formula of Ref. 27J is 
ujq = 10.45). In order to estimate the accuracy of the in- 
tegral over the negative part, let us calculate the integral 
— dujac{ui)/(uj 2 2n 2 a(Za) 2 /2s) so that the terms of 
higher orders in l/u> are accounted for in ac{u>) one af- 
ter another. The results are presented in the table (|TT|) . 
One can see that the successive terms from Eq. (Tj^j) thus 
integrated give the contributions of the same order, i. e. 
the process does not converge to a certain value of the 
integral. 

It is also quite possible that, in order to resolve the 
contradiction between the results ([55]) and ([4"0"j) , the pair 
production in bound-free states should be taken into ac- 
count because of the strong cancellation Eq. (j4"Tj) of the 
contribution of free-free states. 

The expression (j40|) coincides with that calculated 
in Ref. [32j (more precisely, — -D1//2B in the nota- 
tions of Ref. |32j). The comparison of our results, i.e. 
(h(Z) - hs) /hs, and the results of Ref. '32} -D y / f 2B 
is made in Fig. [5] 

It should be noted that our results and those of Ref. [32| 
are essentially different because the last one have a non- 
monotonic dependance on Z. 



j-FACTOR OF A BOUND ELECTRON 

The amplitude of virtual light-by-light scattering is 
known to be a part of the so-called magnetic loop contri- 
bution to the g-factor of a bound electron Ref. [34[ . 

For the lS'1/2 electron state, this contribution reads 
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FIG. 5: Our results (triangles) for (f2(Z) — /2s) //2s and the 
approximation formula Eq. (|31|l (dashed line) in comparison 
with the results of Ref. [32| (squares) 



0.4 



0.2 



FIG. 6: The squares represent the difference Ag Ro f.[2iJ — 
the solid line corresponds to the function 
(16/3) a (Za) 5 (3.35 • 10~ 4 (Za) 2 ), the dashed line is the func- 
tion (16/3) a (Za) 5 (3.35 ■ 10" 4 (Za) 2 + 1.6 ■ 10" 4 (Za) 4 ), cor- 
responding to the Coulomb corrections to the form factor /1 
Eq. (H)}. 



(see Ref. |33(): 



32a(Za) 2 f 00 , r . . . 



j dr r/j (r) f 2 (r) ( ^l^- 



cos qr 



(43) 



where /j and. / 2 are the radial parts of the electron wave 
function in a Coulomb field: 



Jin 

-if 2 (a ■ n) fi 



(44) 



where fi is the spherical spinor (see for example Ref. [35j). 
Using the lowest-order Born approximation for the form 
factor fi and the nonrelativistic expressions for the com- 
ponents of the wave function fi (r) = 2exp(— r/as) and 



f% = /((r)/2m yields the leading correction to the g- 
factor of a bound electron in ISin state Ref. 34j : 



Ag = a(Zaf. 

y 216 v ; 



(45) 



In the case of small Za, one can expand the form factor 
fi in power series of Za: 



/i(0,0,q, Za) = F 



<l 



(Z a ) 2 F {1) (^)+0(Z^). 

(46) 



The contribution of F (q/m) was considered in Ref. [33j 
in detail. To calculate the correction in Za, it is suffi- 
cient to use the functions /1 and /a in the nonrelativistic 
limit and the expression pSD]) A -7/1152 as F(i)(0)- The 
results of the numerical calculation of the magnetic-loop 
contribution exact in Za are presented in Ref. [2l( . 

The comparison of the contribution of the Coulomb 
corrections to the form factor fi Eq. (|30|) and the dif- 
ference of the results obtained in Ref. [21( and Ref. [33| 
is depicted in Fig. [6] It is surprising that our correction 
coincides with this difference not only for the small Za, 
but for Za ~ 1 also. 



CONCLUSIONS 

The Coulomb corrections to the Delbriick scattering 
amplitude have been considered in this article. We have 
calculated these corrections in the low-energy limit but 
taking into account all orders of the parameter Za. The 
accuracy of the calculation is of the order of one percent 
for Z = 50 and increases with Z. Our result is in a 
good agreement with the corresponding contribution to 
the (/-factor of a bound electron calculated previously 
in Ref. [2l|. However, there is a contradiction with the 
dispersion integral of the Coulomb corrections to the pair 
production cross section calculated in Ref. [32| . 

We would like to thank A.I. Milstein and R.N. Lee for 
their helpful comments and discussions. 



EXAMPLE OF REARRANGEMENT OF A 
CONDITIONALLY CONVERGENT INTEGRAL 

While calculating the contribution of the first order in 
Za, we have expanded the expression ([5]) on Za and in- 
tegrated over s. The equation (f2"2"| corresponding to the 
contribution of the diagram [Fig.QJb)] has been derived 
by integration over p, a and t in order [see Eq. (|20|) ]. 
However, one can calculate this contribution in an alter- 
native way - by substituting the expansion of Eq. ^ in 
Eq. (jlip and integrating over r\ and ri before the inte- 
gration over s\ and S2 in the Green functions (0. One 
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of the typical expressions appeared is 



Y(h,t 2 ,z) 



(l + 2z(< 1 + i 2 ) 2 r(t 1 +t 2 ) 
6(1 + M 2 ) S 



x 1 



(h + t 2 y 



, (A.47) 



where t\ 2 = coths 12 G (1, 00) and z = (l + x)/2 S (0, 1). 
The expression (|A.47[) has to be integrated over the total 
variables' domain. One can easily integrate over t\, z 
and t 2 one after another (or z, t± and t 2 ) and obtain a 
finite result, that is 



00 />! foo 



Y dti dz dt 2 : 
1 Jo Ji Ji Ji Jo 



OO fOO r- 1 



Y dz dt\dt 2 



133 13tt 119arctan2 v / 2 38, 32 

= H r- -\ log — . A.48) 

60 4V2 60^ 15 8 81 v ' 



However, if one integrates Eq. (|A.47[) over t\ and t 2 at 
first then the result is the function of z 



00 / />oo 



Y(z) = j (J Y(t u t 2 ,z) dt x ) dt 2 

1 r 1 + i20z 



60z 1 2V27 



^arctan ^2V2z^ — — ^ 



+16z(5-48z)log ( _JL^ ) , ] 
8z 



which has a nonintegrable singularity at z — 



f(z) 



240V2z 3 / 2 



(A.49) 



(A.50) 
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